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XIII. 

A New Method for Correcting a Planets Orbit. 
By TRUMAN HENRY SAFFORD, A.B. 



( Communicated by Professor Peirce.) 



1. It will be assumed that values of the elements are so nearly known, that the 
squares of their deviation from the true ones can be neglected ; the mass of the planet 
will also be neglected. 

2. The approximate elements will be denoted as follows : — 
M° the mean anomaly at the origin of time ; 

p° the semi-parameter of the orbit ; 

tp° the angle of eccentricity (e° = sin cp°) ; 

7r° the longitude of the perihelion ; 

Jl° the longitude of the ascending node ; 

i° the inclination ; 

L° (the mean longitude at the time t a ) = M° -\- n° ; 

a° (the mean distance) = p° sec 2 <p° ; 

k 

[i° (the mean daily motion) = -^— a°~*, where log k == 8.2355814 — 10 ; 

log k — log sin 1" = 3.5500066. 

The corrected elements to be found will be denoted by the same letters without the ° ; 
and the prefix S to any of these letters denotes the correction to be applied to, or any 
variations of, these same elements. 

Analogous modifications of all functions of these quantities will be expressed in the 
same way. 

"We shall use the expressions D M , &c, however, to denote the partial differential co- 
efficients with respect to M, &c. of their functions, even when, for instance, M has the 
value M° ; and so with the rest. 

3. But instead of computing directly the partial differentials of any coordinates with 
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respect to the elements, we shall introduce certain intermediate linear functions of the 
variations 8 M, &c, which last may be considered as themselves variations of conceiv- 
able finite functions of M, &c. 

4. There are, however, two distinct kinds of changes to which the elements are 
subject; we may (1) vary the magnitude of the orbit and its position in its own plane, 
or (2) the position of the plane of the orbit. If the latter change alone takes place, it 
will affect i2 and i, and through £2 the value of tc ; so that in this case 8 it = 8 £2 
(1 — cos i). But n is also subject to other changes, those of the kind (1); it will 
therefore be better if we make 

8 x = 87t— (1 — cos«) 6 si. (l) 

Putting, for convenience, %° = n°, and substituting % = %° -j- 8 %, n = n -j- 8 tc , ■ 

X = n — (1 — cos i) d SI. (2) 

We shall, instead of the longitude n of the perihelion, substitute the distance % of 
the perihelion, counted from that point, A, fixed in the orbit, which is as far back, in a 
direction contrary to the motion, from the ascending node, as the equinox is; the 
former distance being counted upon the orbit and the latter upon the ecliptic. 

5. Let now r, v denote the planet's actual radius-vector and true anomaly, at the 
time t ; and (as before stated, see § 2) r°, v° the same, calculated with the approximate 
elements ; and S r, 8 v the variations r — r°, v — v° ; we shall then get these quantities, 
as also the planet's heliocentric longitude, A, counted from the point A, just alluded 
to, in the following way : — 



E° — e° sin E° = ju° t + M°, 
r° cos v" = a (cos E° — e°), 

r° sin v° = a" cos <p° sin E°, 

jf = v° + n°. 



(3) 



E — e sin E = /t t -j- M, 

r cos v = a (cos E — e), 

r sin v =a cos g> sin E, 

A =zv-\- x . 



(4) 



6. Let now z be a new function of the elements and their variations, and of the time 

t, such that 

€ — e° sin © = ii° z -f M°, 

X cos t) = a" (cos ffi — e"), 

X sin tJ = a° cos <p° sin Q£, 

A = V + n, 



(5) 



We have thus obtained A, the direction in the orbit-plane as seen from the sun, as a 
function of the approximate elements and s, precisely as it is of the corrected ones and 
of t ; but in order to determine the position in space which the varied elements would 
give us, we must superpose other variations. 
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7. The first of these will he 8 w, such that if w° = 0, 

r = tc»-""=: tc iw , (6) 

c heing the hase of the Napierian system of logarithms. 

By its aid, and that of z, we can take account of all the variations of the elements 
which give the dimensions, &c. of the orbit, but not of changes of the orbit-plane 
itself; and z may be considered as the time when the planet would be in the same 
heliocentric longitude (modified as before, like n), as it really is, if its elements were 
M°, &c. The other variation, 8 w, may be denned as the change necessary to adapt 
the natural logarithm of r, calculated by means of the approximate elements and for 
the time z, to its true value for the time t. t 

8. It will be noticed that we have so far adopted a similar course, in representing the 
effect of changes in the elements, to that which Mr. Hansen has employed in his theory 
of perturbations. Our mode of expressing it is taken, with the necessary changes, 
from Prof. Zech's article on Hansen's Method, Vol. XL. of the Astronomische Nachrich- 
ten, reproduced by Prof. Encke, Mathematische Abhandlungen der Berliner Akademie 
for 1855, p. 39. 

9. The orbit-plane will be made to assume a new position by the change of Jl and i ; 
and the planet will be at a distance 8 Z from its position in the unchanged plane ; the 
axis of Z being perpendicular to the orbit. 

10. Our first work will be to compute the partial differentials of z, w, Z with respect 
to the elements. It may be noticed, however, that these quantities, z, w, Z, enjoy but 
a temporary existence as functions of the elements ; they correspond to nothing 
actually existent. We shall, for convenience, omit the °. 

11. By the Theoria Motus, article 15, p. 15, 

7 a a cos op 7 .-_ . . , aa(r-\-p) . _, , 

dv = d (M.4- ut) -\ v ~ rj sin E d q> ; 

r r ' ' rr a 

Hence, d^^^-^^+^sinE^ + rf^ 

Letting now L — L, = it — %, and L, — M = %, there will result 

, , a a cos q> . /T . .\iOor-l-p.—,,,, ,_. 

^■^=-77^^(L,— X + t>t)+— -3=£nnE «** + «**. (7) 

From (5) we get SA = 8z d ^; (8) 



and as z° = f and 8 z is very small, we can assume -j- = -j- . So will 



dz 

mo 

dz dt 



SA=dz d -£. (9) 
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The law of areas gives us 

dA 7 / 

rr— = k\/p = aa cosy./i, (10) 

and from (7), (9), (10), 

3g = 8^ + ^ + Jl+P_ sinE3 + i/_J^_ 1 \, (n) 

ix ' li ' fi a cos q> * ' p \a a cos <p / v ' 

It may be sometimes more convenient to use 

ti8z = d~L, + tdfi+ l^JL s i n E d q> + ('— — l)d x . (12) 

r ' r ' o cos qp ^ ' \a a cos <j> / y 

/ Z* \ 3 

12. The equation p = \j^j,) i« _? cos2 <p 

gives logp = +§ log A — § log sin 1" — f log ,< + 2 log cos <p ; 

and therefore, § log ^ = -j- f log & — f log sin l" — log p + 2 log cos <p ; 

-— = — frflogp — 3 tan <pdcp. (13) 

If, as we propose, we consider /xasa function of log p and 90, we shall get 

Di ogP z = — #* (14) 

D<t z = — 3 1 tan a> + r '~ JJ sin E, (15) 

Y T ' (iflCOSIf v/ 

which will be combined with Dl, « = - ; (16.) 

D % z = - (—^- lY (17) 

z (i \a a cos <j> / ' ^ 

13. As is well known, »• = 1+e J (A _ x) s (18) 
but we also have, by (6), r = r c tt , 

rflog r = <Hogr + ^w, (18') 

, , d log x j 

dlogt = -^-dz, 

by considerations mentioned in § 11. And as, likewise, we can omit the products inter 
se of the variations of r, among which is v — ■ r, we shall have, near enough for our 

present purpose, 

d log r d log r 

~dT~ dt ' 

SO that <5w = dlogr — Slogr, 

-1 * d log r 

= 31ogr — 9x-£-, 

d \o& T 
= dlogr — dA-jf r ; (18") 
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And hence, by (18), (18"), 

„ , cos op cos (A — v) . e sin (A — y) 

8 w = 6 log » — — j-i )-. — *' 8q> — -—. v . . * J . 8 x , 

or l-\-ecos(A — x) l-|-ecos(-4 — x) 

„ , r cos <p cos » „ e r sinv . 

= siog P z — 5 9 — a*; 

J" COS t? 

= 3 log p 8 w — tan m sin E 8 y ; ('19') 

° ^ a cos a? ^ T v ' 

or, =-!^-(2tan<p + ^^)3 9 -tan 9 sinE3*. (19') 

If the planet's place in the corrected orbit be referred to the plane of the uncor- 
rected, it will not be in this plane, but above or below it, if the node and inclination 
do not remain unchanged. The coordinate 8 Z, perpendicular to the plane (A), which, 
as before stated, is that of the approximate orbit, can be ascertained by the usual 
method. 

Thus, let (counted from the ascending node upon the ecliptic) be the ascending 
• node of the corrected orbit upon the plane (A), and 8 I its inclination to that plane. 
We shall then have 8 Z — r sin 8 I sin (v -f- — 0), where a = n — SI. 

But to determine and 8 I we shall have, by spherical trigonometry (omitting 
infinitesimals of the second order), 

sin 8 I sin = sin i sin 8 SI , ) 
sin 8 I cos = sin 8 i, ) 

Or, 8 Z = r sin 8 i sin (v -\- eo) — r sin i sin 8 SI cos (y -j- a>). (20) 

14. So far, we have briefly shown how to resolve the effect upon the planet's place in 
its orbit of slight variations of the first order in the elements, into three portions : the 
first, which, divided by its velocity at the time, we have called 8 s, we have considered 
as the variation of a function, s, of the time t; differing from t only by variations of 
the order of the changes in the elements. This is, of course; in the direction of the 
tangent to the orbit. The second portion we denote by r 8 w, and is in the direction 
of the radius- vector ; the third, perpendicular to the plane of the orbit. 

The first practical use we shall make of this will be to compute the variations of any 
geocentric coordinates with respect to the elements ; and afterwards study the same 
relations with respect to certain new geocentric coordinates which we shall adopt. 

Goetze, Erganzungsheft su den Astronomische Nachrichten, Altona, 1849, has shown 
that if the planet's geocentric place be defined by its position with regard to any 
plane (II.), and if any other plane (I.) (with regard to which o , i2 , i , and % denote 
what the distance of the perihelion from the ascending node upon the ecliptic, the 
longitude of that latter point, the inclination, and, lastly, the angle ^, used above, do 
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with regard to the ecliptic) be employed as a plane of reference for the elements ; and 
again, if I b denote, with respect to the other plane (II.), what geocentric longitude and 
latitude do to the ecliptic ; if, yet again, R, S, W denote, with Encke, the directions 
respectively of the radius-vector, the perpendicular to it in the direction of the motion, 
and the perpendicular to the orbit ; and also F, the direction of the planet as seen from 
the earth ; E, a direction perpendicular to F in a plane perpendicular to (II.), the 
" longitude," I, of any line in which is the same as that of F ; and finally, D, a direc- 
tion perpendicular to both E and F (the positive direction of E being that of b = 90°, 
and that of D, that of ? -|- 90°) ; we shall have (using, with Encke, R S to denote 
the angle made between R and S, and so with other letters), 

A cos b 8 I = cos E D 8 r -f- r cos S D (8 v -\- 3 Xo ) + 8 Z cos W D, 

A 5 b = cos E E 8 r -f- r cos S E (3 v + 8 Xo ) -\- 3 Z cos W E, 

8 A = cos E F 8 r -\- r cos S F (3 v + 8 Xo ) -f- 3 Z cos W F, 

In these formulae we have put (see (20) above) our 8 Z for its equivalent, 

r sin (v -\- a>) sin 8 i — r cos (v -f- <a) sin i 8 SI, 

which is a distance in an absolutely determined direction, and must therefore be the 
same if we employ » , i , I2 , instead of a, i, SL. In general, however, the plane 
(I.) will be the ecliptic. If we transform the above-cited formulse, as Goetze has 
done, into the following, 

cos E D = sin W D sin M' cos E E = sin W E sin N' cos E F = sin W F sin P' ; 
cos S D = sin W D cos M' cos S E = sin W E cos W cos S F = sin "W F cos P' ; 

and make 

D £ r = D £ r = q a cos Q°, 



Erg. Heft, 
(4), p. 162. 



whence, - 



»-I>«^ + ^o) = ? ,> sinQ», 



q" cos Q° = -j- e sin v = — ±£- tan 9 sin E ; 
? «sinQ» = ^ ; 



we shall have, by similar equations to his (11), (12), in the article above cited, 

cos I D z l — 1- sin W D sin (M' + Q°), 

D z b = q ~ sin W E sin (N' + Q°), 
D, A = q" sin W F sin (P' + Q°). 
Again, as d w r = r, 

T> a v = 0, 

the equations before quoted give us, 



(21) 
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And, finally, 



cos b D w I = - sin W D sin M' ; 
A 

D w b = ~ sin W E sin N' ; 
D 8 J = rsin¥F sin P'. 
cos WD 



(22) 



cos b Dz I 
B z b = 



A 

cos WE 



(23) 



Dz A — cos W P. 

The following are (Astr. Nachr. XXVIII. 115) Goetze's general expressions, by 
means of which W D, M', &c. may be computed : 

sin WD6inM = — sin (/ — Jl') ; 
sin W D cos M = cos (I — Jl') cos I ; 
cos W D = — cos (I — Jl') sin I ; 
M' = M + a)' + i7; 



(24) 



sin W E sin N = — cos (I — Jl') sin b ; 
sin W E cos N = — sin (I — Jl') sin 5 cos I -f- cos b sin I ; 
cos W E = sin (I — Jl') sin b sin I -\- cos b cos I ; 

N' = N + w' + v ; 



(25) 



(26) 



sin W F sin P = cos (I — Jl') cos b ; 
sin W F cos P = sin (I — Jl') cos b cos I -f- sin 5 sin I ; 
cos W F = — sin (I — Jl') cos b sin I -\- sin b cos I ; 
P' = P -f <»' + v. 

In these expressions, £2', o, I, denote respectively the distance of the planet's as- 
cending node upon the plane II., counted from the same origin as I is, the distance (in 
the direction of motion) from this ascending node to the perihelion, and the inclination 
to the plane II. 

15. We will now let the approximate orbit-plane itself become the plane II. We 
have first to devise means for referring observed geocentric places to it. We shall, 
therefore, suppose that the planet's geocentric orbit-longitude, so to speak, is denoted 
by r\, it being counted from the ascending node upon the equator ; that 6 denotes, like- 
wise, its geocentric orbit-latitude, referred in the same way to the orbit. 

It must be observed that throughout this portion of the investigation we are speak- 
ing of the approximate orbit as a fixed plane of reference. 

Let Sl x , Oj , ii , denote the approximate longitude of the ascending node upon the equator, 
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the distance (in the direction of the motion) of the perihelion from this node, and the in- 
clination to the equator. Let a, S denote the geocentric right-ascension and declination. 

Then will tf, 6 be the same functions of a ■»— J2 X , d, i u that the geocentric longitude 
and latitude are of a, S, and the obliquity. 

We shall thus have (Theoria Motus, Art. 68, p. 64), 

sin (45° — £ 6) sin £ (E — v ) = cos [45° -f £ (« — .ft,)] sin [45° — £ (*, -f 5)] 
sin (45° — J 6) cos £ (E — v ) = sin [45° -{-£(«_ .ft,)] cos [45° — J (*, — 3)] 
cos (45° — £ 6) sin J (E + n ) = sin [45° + £ (a — J2,)] sin [45° — J (i, — 5)] 
cos (45° — £ 6) cos £ (E + v ) = cos [45° + £ (« — .ft,)] cos [45° — J (t, + 3)] 

It will not be necessary in general to compute -q, 6 with extreme accuracy, as what 
we most wish are their variations. These will be (see same article in Theoria Mdtus, 
adfinem, where, as in (27), E is not the eccentric anomaly), 

d r\ cos 6 = sin E d a cos 8 -j- cos E d d ; 



(27) 



(250 ";, ;;„", "' , \ -(26) 

■7]. ) 



(270 
d 6 = — cos E d a cos d + sin E d S. ) 

In applying now the formulae (24) to (26), we shall have 1 = 0. J2,', indefinite, 
but Jl' -\- ca = el , as the origin of q is at the ascending node, on the equator. The 
letters r\ and are substituted, as the things which they denote are here the same, for 
I and b, Goetze's notation. 

Then, WD = 90° ; . 

M = — t, -f- Jl' ; ( (24') 

M' = w, -\- v — t} ; ) 

WE=6; ) WF = 90° — 6; 

N' = 270° + (»! + » — V ;l ( ' P' = 90° + *>, + *>• 

From (23), (24') we have the expression 

Dz n = ; (27) 

and as the variations g, w do not contain Ji or i, and are thus (to terms of the first 
order with respect to the changes of the elements) independent of small changes in 
the orbit-plane, so also is ij. 

That is, so soon as we know the positions of the orbit-plane correct within limits 
of error of the first order, we can obtain one geocentric coordinate which depends only 
upon four elements. 

16. If a set of rectangular coordinates be employed, the axes of x, y in which are in 
the approximate orbit-plane, and the former passes through the nodes upon the equator, 
the positive direction of x being towards the ascending node, and if S, T, *P denote 
the sun's geocentric coordinates, we shall have, if the elements are supposed unvaried 
(J denoting, as usual, geocentric distance), 
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r° cos (<»i° + i>°) -\-S—A° cos r? cos ff> ; \ 

r° sin (a,, -f v") + T = A° sin r? sin (9° ; \ (28) 

W = A° sin 0°. ) 

The elements being supposed changed, we shall have 8 (g>! -|- v) = 8 A, provided we 
consider 8 ai as the variation arising from a change in the position of the perihelion, 
and do not permit it to be affected by the change in the orbit-plane. 

So that Q being any angle whatever, and &, T, ¥ being invariable, we shall have 

8 [r sin («, + » + Q)] = 8 [A sin ( v + Q) cos 8~\. (28') 

17. We will now apply the formulae we have obtained to a practical problem of very 
frequent occurrence, — the case in which a set of elements is desired which shall satisfy 
two observed or normal positions precisely, and certain others nearly. 

In our case we will, with regard to those positions which are not to be exactly 
represented, content ourselves with satisfying the observed value of our coordinate, rj ; 
the plane to which it refers being that of an approximate orbit satisfying the first- 
mentioned positions. 

We will, therefore, assume values of J corresponding to times t t , &>, for which we 
have observed right-ascensions. From this the orbit is easily found, corresponding to 
these distances. 

A variation of the distance for the time t, which we may denote by J, will produce 
the following change in equation (28) : 

8 [r sin (w, + v -f- Q)] = d A . sin (ij -f- Q) cos 6. (29) 

But we put 8 h t = 8 (ax -\- v) under the form 8 s . D t v = 8 z . D e A, as in (9) ; 
whence by (18) 

S log r = D t log r . 8 z -\- r 8 w, 

and (29) becomes 

8 [r sin (o>i -\- v -\- Q)] = D { [r sin (m, -|- *> + Q)] 8 z + r sin (w, -|- v -{■ Q) 8 w . 

Making now D t r = -~ tan <p sin E = c cos y, } 

> (30) 

r D ( v = — — = c sm w, l 

r ~ J 

we shall have, if 8 denote the effect of changing J only, 

8 A sin (jj 4~ Q) c° s 6 = c sin (ta, -)- v -\- y -f- Q) 8 z -f- r sin (a)! -f- v -(- Q) 8 w. (31) 

If Q = — e^ — v — t//, (31) becomes 

8 A sin ()j — Wi — v — if) cos 6 = — r sin y/ . 8 w ; 

_ sin (o>i + » + ^ — l) /, , M > 

D 4 w= — v ' . ', -* "cos<9. (32) 

r sin i/r ' 
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So also, making Q = — a t — v, we find 

_ sin (n — o>i — v) . r sin (n — on — v) . ,„„. 

Da z = — !l/ — ; — I cos d = -,—r J y cos 0. (33) 

cuay « VT 

If ^/, and consequently the orbit, is changed, so that the plane to which the coor- 
dinates are referred is no longer the actual plane of the orbit, the last equation of the 
group (27) becomes 

S Z -j- v = A sin 6. (34) 

And thus D A Z = sin 6. 

An arbitrary variation in either assumed geocentric distance is thus separated into 
three components ; one, c 8 z, in the direction of the tangent to the orbit, at the planet's 
place for the time to which the varied distance refers, c being the tangential velocity ; 
another, r 8 w, in the direction of the radius vector, and so not making an angle = 90°, 
but one = f, with the former; and a third 8 Z, perpendicular to the orbit-plane, which 
is the plane in which the other two lie. 

A new orbit computed with the distance J t , varied by the amount 8 z/ 1} and the dis- 
tance 4% unvaried, but still using a l5 a 2 , 8 L , 8 Q , will be connected with the former in the 
following way. 

If f x = D Al z t , g x = Da, w t , then there will be four equations, of which the known 
terms will be f x , D A , z 2 = 0, g x , D A , w a = 0, and the unknown the differential coeffi- 
cients of these elements which lie in the orbit-plane ; the coefficients of these unknown 
quantities being the coefficients of the variations of the elements in (11) and (19'). 

JBut instead of fi as an element, we shall use p, whereby 

d log fi = — | d log p -\- 3 d log cos <jd ; or, 

d log ji = — f d logp — 3 tan q> d <p. (13) 

Our four equations will be, then, 

/l = J D . i L / ^^ 1 D ai logp + ( / ^sinE 1 --3tan r ./ I )D Ai < P + i(^-l)D. 1 ^ 

•0 = iD Al L # -S«,D Al ]ogj, + (^^BinB»-8tan 9 .4)D Al9 + i(^-l)D 4| ,; 

M ^ ^ x v ' ) (35) 

•r^ . »"l COS Vl _ . _ _ 

o, = Da, log p Da, a> — tan a> sin ^ Da x ; 

J" I D -' a cos qp l • 

_ . r« cos »j _. , • -n T-. 

= Da log » Da„ <j> — tan m sin E 2 Da, v ■ 

i DJ a cos cp * T T ! 

By subtracting the first of these equations from the second, we can at once reduce 
our equations to the number of three, with three unknown quantities ; determining 
Da L, afterwards from the second equation explicitly, or, what is perhaps better, leav- 
ing - • D Al L, expressed in terms of the variations D A , log p , &c, in which form it can 
be easily used. 
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If now we vary z/>, we shall have similar equations, in which, after performing the 
subtraction indicated above, we shall have, in place respectively of — fi,(ji,0, as 
known terms, the values^, 0,# 2 ; an d a^ unknown quantities, the first derivatives- of 
the elements with respect to A^. 

18. Having now these derivatives with respect both to J x and 4%, we can for any 
other time easily compute, by formula? (11), (19), using the former as modified by (13), 

32 = ^L / -^51ogp + ( / ^+-^sinE-3e.ta ny )^ + i(^-l)3 / , ) (36) 

the values of D Al z, D A , w ; and as, by (26'), 

the value.(which we do not know) of 8 Z will not affect y ; and the coefficients D A[ , 
D^ are, with regard to almost all planets, of such slight amount that they cannot 
increase the weight of our solution much over T ^ 7 part, as has been shown in Mr. G. 
P. Bond's Memoir on " Equivalent Factors," assuming the modulus of D A , as about 
^, which is as much as it can often be with regard to five sixths of the asteroids. 

19. We have thus, if at least two positions besides those to be exactly satisfied are 
given, the means of determining what changes the assumed J lt J z must undergo to 
satisfy all the geocentric places we have, with almost if not quite the accuracy of a 
least-square solution. After this is done, it only remains to compute from the variations 
of Ji, Jz the changes which the orbit-plane must undergo, in order to pass exactly 
through the positions to which they refer. This is done by combining (20) and (34). 

We have 

8 Zj = sin #! 8 z/i = r, sin (y x -J- ui) . 8 i — r x sin i cos (v, 4- to) . 3 SI ; ~i 

> (37) 

d Z 2 = sin 2 d A% = r s sin {v % -j- «) . 8 i — r 2 sin i cos (v 2 -\- w) . 8 SI. ) 

The readiest practical way to obtain d i, 8 £2 will be to form these equations and 
solve them numerically. 

After this is done, we shall find also, by the equations in which D Al , L , &c. occurred, 

the values of 8 L , 8 log p, 8 <p, 8 %; and the resulting changes in the usual elements 

will be, 

8 L = 8 L, -\- (1 — cos i) 8 Si ; 

8 a = 8 x + (1 — cos i) 8 Si ; 

8 log a = 8 logp 4" ^ tan qp . 8 <p sin 1"; 

8 (i = — g (i 8 fogp — 3 fi tan q>.8q> sin 1". 

It is to be noticed that all the variations, those of A x , J%, and log p inclusive, will 
be expressed in seconds; so that we must add to their logarithms the log 4.6855749 — 
10, or must divide them by 206264.8 to reduce them to the unit required. 



(38) 
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20. As an example, we have selected the planet Aglaja @. But few observations 
of it were made near its first observed opposition, in 1857 ; and five normals contain as 
much of the material at hand, as, on account 'of its scattered nature, we have thought 
it well to use. For several very exact normals do not help the matter much, if they 
must, from the nature of the case, be so combined with less exact ones that the burden 
shall fall on these last ; the chain is not stronger than its weakest link. 

The normals are : — 

d 

O i a 

— 49 44.2 

—2 11 17.4 

—1 45 47.6 

+3 32 41.7 

+7 52 1.2 

Assuming the following for the logarithms of two geocentric distances, 

Sept. 17.5 0.201821, Nov. 16.0 0.307251, 

the resulting elements, referred to the equator and mean equinox of 1857, Jan. 0, were, 

Epoch, 1857, Nov. 16.0 



lumber. 
I. 


Date. 
Washington M. T. 

1857, Sept. 17.5 


No. Obs. 
13 


a 

O 1 ft 

21 29.8 


n. 


" Oct. 21.0 


11 


354 21 7.6 


hi. 


« Nov. 16.0 


9 


353 38 16.4 


IV. 


1858, Jan. 11.0 


3 


3 52 10.9 


v. 


" Feb. 10.5 


2 


13 21 32.7 





M° 


46 20 


24.3 






OJl 


313 31 


56.7 






fli 


49 


9.5 






h 


28 27 


7.0 






O 


7 21 


56.8 






o 


12 


5.719 




The other normals were represented 


thus, 






Number. 

n. 


Sate. 

Washington 

1857, Oct. 


M.T. 

21.0 


C- 

Aa 

+ 0.4 


-O. 

Ji 

II 

+ 0.1 


IV. 


1858, Jan. 


11.0 


+25.4 


+11.9 


v. 


" Feb. 


10.5 


+54.9 


+33.0 



For those dates the equations (27), (27') give us, 

r, 6 E • C — O. 

II. 353 15'.58 +1 9.04 61 44.36 d ij = + o!40 5 6 = — o'!lO 

IV. 4 21.78 +1 40.38 61 34.72 +27.96 —1.60 

V. 14 44.20 +1 1.88 62 16.84 +63.50 +3.91 

For the dates I. and III., we find 

I. 359 11.38 —0° 30.50 

m. 352 49.85 +1 51.75 
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As we do not wish E here, we can apply a method like the simpler well-known forms 
of converting right-ascension and declination into latitude and longitude. 
From the elements 

log (k V>°) = 8.461728 ; 

and from what has been previously obtained, 

log r v ip 

Sept. 17.5 0.413888 43° 50.67 85° 21.16 

Nov. 16.0 0.423901 58 11.80 84 10.36 

From (32) and (33) are now derived, 

log/. log g. 

For *, = Sept. 17.5 +0.45230 +9.58475 ) 

< 2 = Nov. 16.0 —1.47254 +9.56672 \ 

The equations (36) for determining D Al , log p, &c. become, subtracting the second 
of them from the first, and also those with the same coefficients for finding D A2 log p, 
&c. if we make D Al log p = l x , D Aj \ogp = ? 2 , D Al q> = m x , D Aa <p = m^, D Al % =■ «,, 

(1.) (2.) 

+2.833 +29.685 = (1.95061) I— (1.83514) m — (1.03991) n; 1 

+0.38437 0.00000 = (0.00000) I— (9.81607) m — (8.91004) n; I (39) 

0.00000 + 0.36874 = (0.00000) I — (9.68983) m — (9.00885) n ; J 

of which the solutions are, 

log (D Aj log nat p) = +9.8222 ; log (D A , cp) = —9.9737 ; log (D Aj X ) = +1.0423 ; 

log (Da 2 log rmtp)= — 9.4000; log (Da 2 9 ) == +9.9620; log (Da 2 x) = —1.0199. 

When J x is changed, the value of z for % remains unchanged, and vice versa ; so 

that, for any time, 

D Aj z = D Aj (2 — 2 a ) ; (40) 

But z — z 2 does not contain L , so that we shall have 

D Al z = Diogp (2 — z 2 ) D Aj cos p + D^ (2 — 2 2 ) D Aj <jp + D x (z — z 8 ) D A *. (40') 

Similarly, D*. * = Da, (* — *,). (41) 

Nor does 8 w contain 8 L,. 

We now find for the times of normals II., IV., and V., by formulae (11), (19'), as 
modified in (40), (41), 

D Al 2 D A( w Da 2 z D A]J w 

H. — 4.373 +.1645 — 11.33 +.2108 

IV. +37.39 —.3285 — 96.07 +6848 

.V. +72.52 —.4890 —146.44 +8395 
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t 

And we thus finally form the equations in which 8tj cos is the known term, 8 J x , 
8 J% the unknown quantities, and cos 6 D Al rj , cos D Aa rf the coefficients. 
From normals 



H. = +0.40 + 0.0249 8 4, — 0.0036 3 A % 
IV. 27.96 0.0284 —0.1259 

V. 63.50 0.1145 —0.2672 



(41) 



From these equations were derived, by least squares, the values 

dj l = — 6 sin l", 

3 4, = +231. 2 sin 1", 

the weights being made proportional to the number of observations. 

By means of the previously computed partial differential coefficients, we are now 
enabled to obtain the values 



' > in the sixth place of decimals. 
',p = — 130.6 > 



3 log nat p — —300.9 
.\ 3 com log j 

logp — 0.452163, 
8 <p — 3' 37".4, 

9 = 7° 25' 34".2, hence, log cos q> 9.996342. 
3 x = x — n" = — 41' 26".6, 
log a = 0.459479. 

Referring the previous, and still unchanged, node and inclination to the mean eclip- 
tic and equinox of 1857.0, we find 

Jl° 4 28 34.7, 

6>° 309 47 34.0, 

i° 5 25.6. 

Hence, x = 313 34 42.1. 

From the value of D Aj z 2 =fa, we obtain for 1857, Nov. 16.0, 

3 h =/ 2 ,8Ji = —0*03328. 

1c iff) 

Consequently, • 8v 2 = 8 J 2 — 3x = ^ sin p 3 « 2 — 3 x- 

This then gives us d » s = — 28".2 — 3 x = +40' 58".4. 

Then Vi = 58° 52' 46".2, 

M = 46° 49' 58".7. 

Check, r t = — p-£ = r° 2 (14-3 wX 

1 + e cos i>2 » v i v> 

which comes out exactly 0.424080 in both. 

From the equations (38) the values of 8 i = -\-l\".% 8 £1 — — 16".l, were rather 
roughly estimated ; whence it — % = 2 sin 2 J i . 8 Jl = — 0".l. 
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The elem 


ients will then be, 






Epoch, 1857, Nov, 


, 16.0. — M. Equinox, 1857.0, 




M 


46 49 58.7 




n 


313 34 42.0 




SI 


4 28 18.6 




9 


7 25 34.2 




i 


5 36.8 




/* 


12 5.752 




log a 


0.459479. 



The remaining errors, S r\ cos 0, after the solution of the conditional equations (41), 
were, 

C. — 0. 

n. — o".6 ] 

IV. — 1".3 l (42) 

V. +1".0 J 

Substituting the new elements in the formulae for the geocentric places, we obtain 
values of 3 a, S 8 corresponding with them ; and from these, by (27), values of 8q cos 6, 
differing from those in (42) very slightly. 

The resulting differences in the other geocentric coordinate, 0, are nearly as follows : 

S6(Q. — O.) 

n. +i".s 
rv. — 2".o 

V. +5".7; 

which is a small difference, seeing that the last normal depends upon but two observa- 
tions, and those, as is apparent by Mr. Oeltzen's investigation (Astr. Nachr., No. 1167) 
disagreeing largely, as indeed do the observations upon which IV. is based. 

Ten months after the date of the last of these normals, the comparison of an ephem- 
eris computed from the above elements with a Cambridge observation, kindly furnished 
me by Prof. W. C. Bond, was as follows : 

C — 0. 

J a cos d A d 

1858, Dec. 28 +1'.6 +1'.2; 

showing an agreement much more satisfactory than is usual in determinations from the 
first apparition of an asteroid, especially considering the weakness of normals IV. and 
V., upon which, as may be seen by (41), almost the entire weight of the computation 
rests. Perturbations, too, were entirely neglected. 

The calculation of the above orbit was executed entirely with six decimal places (less 
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in portions). The usual method (of varying two mean distances, actually computing 
three orbits) would not only have required more operations, but these must, many of 
them, have been carried to seven decimals, in order to obtain with any degree of accu- 
racy the coefficients of the unknown quantities in the equations (41). At least one 
third the labor appears to be saved by this process, while a differential method, in gen- 
eral much better than the rule of false position, is substituted for the latter. 

Moreover, the general relations developed in the first part of this paper are not with- 
out very interesting aspects with regard to the general problem of representing long 
series of geocentric observations by an orbit, by the method of least squares ; and it 
is not impossible that further calculations may be made upon the theory of Aglaja, 
using the numbers and formulae already obtained, and approximating more and more 
closely to the true elements, when the observations requisite shall have been made. 



